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ABSTRACT 

We investigate the prospects for detecting dust from two body collisions dur¬ 
ing the late stages of planet formation at 1-150 AU. We develop an analytic 
model to describe the formation of a dusty cloud of debris and use numerical 
coagulation and n-body calculations to predict observable signals from these 
events. In a minimum mass solar nebula, collisions of 100-1000 km objects at 
distances of 3-5 AU or less from the parent star are observable at mid-infrared 
wavelengths as bright clumps or rings of dust. At 24 /mi, the clumps are ~ 0.1-1 
mag brighter than emission from dust in the background debris disk. In edge-on 
systems, dusty clumps produce eclipses with depths of < 1.0 mag that last for ~ 

100 orbital periods. Large-scale surveys for transits from exosolar planets, such 
as Kepler, can plausibly detect these eclipses and provide important constraints 
on the terrestrial environment for ages of < 100-300 Myr. 


Subject headings: planetary systems - solar system: formation - stars: formation 
- circumstcllar matter 


1. INTRODUCTION 

Throughout history, collisions have shaped the Earth. Repeated mergers of 10 km and 
larger objects produced the proto-Earth in 10-30 Myr (Wetherill & Stewart 1993; Ag- 
nor, Canup, & Levison 1999; Chambers 2001; Yin et al. 2002). A few Myr later, a giant 
collision with a Mars-sized planet removed the Erst terrestrial atmosphere and led to the 


formation of the Moon (Hartmann & Davis 1975; Cameron & Ward 1976; Benz, Slat¬ 
tery, & Cameron 1986; Canup 2004a,b). During the Late Heavy Bombardment rs_/ 4 Gyr 
ago, a sustained period of large collisions shaped the surfaces of all terrestrial objects (e.g., 
Tera, Papanastassiou, & Wasserburg 1974; Hartmann 1980; Neukum 1983; Ryder 2002; 
Koeberl 2003). Giant collisions also produced structure in the Zodiacal dust band and 
perhaps mass extinctions on the Earth (e.g. Alvarez et al. 1980; Nesvorny et al. 2002b). 

Understanding the formation and early evolution of our solar system requires unraveling 
this collisional history from current images, orbits, and samples of asteroids, meteorites, 
planets, and satellites. Although these data yield a remarkably detailed picture of the 
formation of the solar system, the record is incomplete, particularly at the earliest epochs 
(e.g., Swindle 1993; Melosh, Vickery, & Tonks 1993; Wadhwa & Russell 2000). Some of 
the early collisional history is also erased by subsequent events, including volcanic activity 
and collisions (e.g., Wilhelms 1987; Strom, Schaber, & Dawsow 1994). The ambiguity of 
the early record leads to uncertainties in our understanding of the formation of chondrites 
(e.g., Itoh & Yurimoto 2003; Wood 2004; Connolly 2005), the origin of the Late Heavy 
Bombardment (e.g. Cohen, Swindle, & Kring 2000; Levison et al. 2001; Koeberl 2003), 
and other aspects of the early history of the solar system. 

Observations of debris disks surrounding other stars provide another way to study col¬ 
lisions during the evolution of a solar system. At least 40% of all A-type stars and 10% 
of all solar-type stars have far-infrared excesses from circumstcllar dust (e.g., Backman & 
Paresce 1993). Debris disks with masses of ~ 0.01 M® in small, 1-100 fin i grains explain 
the scattered-light images and the spectral energy distributions of the best-studied systems 
(e.g., Backman & Paresce 1993; Artymowicz 1997; Greaves 2005, and references therein). 
However, collisions and radiation forces remove small grains on timescales, < 1 Myr, much 
shorter than the 10 Myr to 1 Gyr age of the parent star (Aumann et al. 1984; Backman & 
Paresce 1993; Zuckerman 2001; Thebault, Augereau, & Beust 2003). If the mass in 1-10 
km bodies is ~ 10-100 M®, debris produced from high velocity collisions among these larger 
bodies can replenish the small grain population for 100 Myr or longer (Habing et al. 2001; 
Kenyon & Bromley 2002; Dominik & Decin 2003). Detailed numerical calculations suggest 
that the formation of debris disks coincides with the growth of 1000 km or larger planets 
throughout the disk on similar 10-100 Myr timescales (e.g. Kenyon & Bromley 2004a,b). 
Thus, debris disks also provide a signpost for the formation of solar systems. 

Recent theoretical work helps to link our understanding of the collisions in debris disks 
with the collisions that shaped our solar system. For debris disks around A-type stars, 
Dominik & Decin (2003) and Kenyon & Bromley (2002, 2004b) showed that the ratio of 
the disk luminosity L c i to the stellar luminosity L* should follow a simple power-law relation, 
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Lrf/L* « L 0 /(l + t/to), where L 0 ~ 10 -3 to 10~ 2 and to ~ 100-300 Myr. Observations 
of infrared excesses in nearby stars follow a similar relation (e.g. Kalas 1998; Rieke et al. 
2005; Greaves et al. 2005). These results suggest an exponential decline in the collision 
rate, similar perhaps to the monotonic decline in the background accretion rate following 
the formation of the Moon 1 (e.g., Ryder 2002; Koeberl 2003). Among others, Wyatt & 
Dent (2002) noted that large collisions are frequent enough to produce dusty clumps in 
debris disks, with long-term structure similar to the local Zodiac (Nesvorny et al. 2002b). 
In the terrestrial zone, Kenyon & Bromley (2004a) demonstrated that collisions between 
50-100 km and larger objects can raise the debris disk luminosity above the simple power- 
law relation, suggesting that moderate and giant impacts are observable for short periods of 
time (see also Stern 1994; Grogan et al. 2001; Nesvorny et al. 2002a; Zhang & Sigurdsson 
2003). 

Here, we continue to develop the idea that giant collisions produce observable structures 
in debris disks. We derive an analytic model to identify collisions that produce observable 
clumps of dust and we use numerical calculations to explore the conditions where large 
collisions are visible. In a minimum mass solar nebula, collisions of 100-1000 km objects at 
distances of 3-5 AU or less from the parent star are observable as bright clumps or rings of 
dust. In disks observed close to edge-on, the optical depth is often large enough to produce 
eclipses of the parent star. These eclipses last for several days and recur over 10-100 orbital 
periods. We use detailed IV-body simulations to illustrate the evolution of these eclipses 
and to make estimates for the frequency of these events. Current and planned surveys for 
planetary transits can detect these eclipses and test our predictions. 

We develop the analytic model in §2, describe numerical simulations of an ensemble of 
colliding planetesimals in §3, and derive the time evolution of single catastrophic collisions 
in §4. We conclude with a brief discussion and summary in §5. 


1 We distinguish a general monotonic decline in the impact rate 4-4.5 Gyr ago from the Late Heavy 
Bombardment, an apparent spike in the impact rate ~ 600 Myr after the formation of the Earth and the 
Moon (see, for example, Hartmann 1975, 1980; Ryder 1988; Hartmann et al. 2000; Ryder 2002; Koeberl 
2003, and references therein). 
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2. ANALYTIC MODEL 
2.1. Derivation 


We begin with an analytic model for the formation of a clump of dusty debris from two 
colliding planetesimals. Each planetesimal has mass density p. The larger planetesimal is the 
‘target’ and has mass rn t , radius r tl and orbital eccentricity e t and inclination i t . The smaller 
‘projectile’ has mass m P) radius r p , eccentricity e p , and inclination i p . The planetesimals lie 
in a disk with surface density E(a), where a is the distance from a star with mass M* and 
luminosity L*. 


We assume that the collision produces a spherical cloud of debris with radius r c and 
mass ■ rrid■ The debris consists of particles with sizes ranging from 1 pm to 1 m or larger. We 
follow previous investigators and derive the debris mass from the ratio of the impact energy 
Qi to the disruption energy Q d , 


rn d = 0.5 (m p + m t ) 



( 1 ) 


where (3 e is a constant of order unity (see Davis et al. 1985; Wetherill & Stewart 1993; 
Kenyon & Luu 1998; Benz & Asphaug 1999, and references therein). If V is the relative 
velocity of the two planetesimals and V e is the escape velocity of a single body with mass, 
m = nip + m t , the center-of-mass collision energy is (Wetherill & Stewart 1993) 


_ m p m t Vj 
1 4 m 2 


( 2 ) 


where the impact velocity Vf = V 2 + V) 2 . The energy needed to remove half of the combined 
mass of two colliding planetesimals is 

/ r \Pb / r \Pb 

Qd = Qb(~ -) +pQ 9 (- -] , (3) 

VI cm/ VI cm/ 

where is the bulk (tensile) component of the binding energy and pQ g r is the gravity 

component of the binding energy (see, for example, Davis et al. 1985; Wetherill & Stewart 
1993; Holsapple 1994; Benz & Asphaug 1999; Housen & Holsapple 1999). 


To derive the amount of stellar radiation intercepted by the debris cloud, we define the 
radial optical depth through the cloud as 


r 


3 f n(r)r 2 dr 


(4) 


where n(r ) is the number of particles with radius r within the debris cloud. The radiation 
intercepted by the debris is roughly « rr 2 /4a 2 . As orbital shear broadens the debris 
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cloud into a ring, r drops but Z^/L* remains roughly constant. If n(r) = no'r~ a and rn d = 
(47rp/3) f n(r)r 3 dr, L d / L* is a function of a, m d , and a and does not depend on the initial 
radius of the debris cloud r c . 


The collision rate places strict limits on a. For the late stages of planet formation, the 
timescale for a collision between an object and smaller bodies in the disk is (Lissauer 1987; 
Wetherill & Stewart 1993) 


t_c 

P 


pr 

7m E 


(5) 


where P is the orbital period. If we integrate the expressions for r and rn d over r, adopt p 
= 1-3 g cm -3 , and set the surface density of small particles in the background debris disk 
to E = E{,a -3 / 2 with E& ~ 0.1 g cm -2 at 1 AU, we can rearrange equation (5) as 


tc 

P 


r^j 


10- 25 r 3 


1 AU 
a 


1/2 



( 6 ) 


In each orbit, particles with t c /P <0.1 have many collisions with material in the background 
debris disk. Because the background has a wide range of sizes, r ~ 1-10 pm to r ~ 1 km, 
objects in the debris cloud rapidly reach an equilibrium size distribution with a ~ 3.5 
(Dohnanyi 1969; Durda & Dermott 2002; Williams & Wetherill 1994). For an observable 
debris cloud with L d /L± > 10~ 5 , particles with r < 10-20 km reach this limit at a = 1 
AU. Although larger particles may have a steeper size distribution (Campo Bagatin & Petit 
2001; Michel et al 2001, 2002), these contribute little to the optical depth. Thus, we assume 
a ~ 3.5 for all particles in the debris. 


With a known, the optical depth and luminosity are 


m d 

/y _ 

np(r)r 2 c 

where (r) is a typical radius of a particle, and 


L d rn d 

_ _ 

L* 47t p(r)a? 


(7) 

( 8 ) 


The reprocessed luminosity is independent of the geometry of the debris cloud and depends 
only on the location of the collision and the mass distribution of the particles in the debris. 
For L d /L$ > 10~ 5 , this expression constrains the total debris mass: 


m d > 10 


21 


(r) 


1 mm J \ 1AU 


g- 


(9) 


Thus, an observable debris cloud requires the complete disruption of (i) two 100 km objects 
at 1 AU, (ii) two 300 km objects at 10 AU, or (iii) two 1000 km objects at 100 AU. From 
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equation (5) these debris clouds form every ~ 10 2 yr at ~ 1 AU and every ~ 1-3 x 10 6 yr 
at 30 AU. 

The timescale for orbital shear to broaden the debris cloud into a ring and mix with the 
background disk is t m « ( a/r c )P. After the collision, the debris has a velocity dispersion 2 
Vd/V k ~ lO -3 (e.g. Michel et al 2001, 2002), which implies an approximate size r c /a ~ 
Vd/V k ~ 1CT 3 . Gravity sets a lower limit to the extent of the debris cloud. Debris cannot 
remain bound to the largest objects beyond the mutual Hill radius, Rn/a = (2m p /3M*) 1 / 3 
~ 1 — 4 x 1CT 4 for 100-300 km projectiles. The disk scale height, H « 0.01-0.la provides a 
rough upper limit to the width of the debris cloud. These estimates suggest that the cloud 
remains distinct for at least 10 2 orbital periods. 

Our approach ignores the increase in thermal energy during a collision. From detailed 
numerical calculations (e.g., Canup 2004b), large impacts produce considerable increases in 
the temperature of the colliding bodies. The results suggest that the temperature of two 
large bodies rises to ~ 5000 K to 10,000 K during the impact (see also Stern 1994). At these 
temperatures, the timescale to radiate the impact energy is short, ~ days, compared to the 
lifetime of the debris cloud, ~ years. If the increase in thermal energy is large compared 
to the impact energy, our analysis overestimates the amount of debris and therefore yields 
upper limits to the debris cloud luminosity. 


2.2. Application 

To apply the analytic model to debris disks, we use parameters appropriate for the late 
stages of planet formation. We adopt a minimum mass solar nebula model for the surface 
density, E = E 0 (a/ao)~ 3,/2 , where ao = 1 AU and E 0 = 15 g cm -2 (Weidenschilling 1977; 
Hayashi 1981). For this section, we ignore the change in surface density at the ‘snow line,’ 
where the surface density increases due to ice condensation. This approximation decreases 
our collision rates by a factor of two, but has no impact on the luminosity or optical depth 
of a debris cloud. 


2 The range in the velocity dispersion is a factor of ~ 2-3. After the collision, dynamical interactions with 
other objects and radiative processes such as the Yarkovsky effect can modify the velocity dispersion (e.g., 
Michel et al 2002, 2004). 
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The size distribution of projectiles and targets is (Kenyon & Bromley 2004c) 


n(r) 


ns'r~ as r < rq 

< nj r i <r< r 0 

n L r~ aL r > ro 


( 10 ) 


Numerical simulations of planet formation suggest r 0 ~ 1-100 km, with ck^ ps 4 for large 
objects and as ~ 3.5 for small objects (e.g. Stern & Colwell 1997; Kenyon & Bromley 
2004a). We use r 0 ~ 10-100 km and rq « 0.1-1 km (Kenyon & Bromley 2004c). For the 
orbital parameters, we use M* = 3 M 0 , e = 0.2, and i/e = 0.5. 


We assume standard values for the bulk parameters of large objects. In our model, 
rocky objects have p = 3.0 g cm -3 , Qb = 6 x 10 7 erg g _1 , ab = —0.37, Q g = 0.4 erg cm -3 , 
and a g = 1.36; icy objects have p = 1.5 cm -3 , Qb = 10 3 erg g _1 , atb = 0.0, Q g = 1.4 erg cm -3 , 
and a. g = 1.25. For this paper, we ignore a weak relation between the disruption energy and 
the impact velocity (Housen & Holsapple 1990, 1999; Benz & Asphaug 1999). 


To explore the sensitivity of our results to the bulk properties, we consider two alter¬ 
natives. For rocky objects, our simple approximation to the empirical Qd derived from the 
size distribution of asteroids in the solar system (Durda, Greenberg, & Jedicke 1998) is Qb 
= 5 x 10 4 erg g _1 , aj, = 0.0, Q g = 1.7 x 10~ 8 erg cm -3 , and a g = 2.4. As an approximation 
for icy rubble piles, we also consider Qb = 1 erg g _1 , ab = 0.0, Q g = 10“ 3 erg cm -3 , and a g 
= 1.25 (e.g., Lcinhardt & Richardson 2000, 2002). 

To calculate the properties of debris clouds, we use the orbital parameters and bulk 
properties to derive rn c i from equations (l)-(3). For n(r) = n 0 r~ a and a = 3.5, n 0 = 
3rrid/(4:7ip f r 3 ~ a dr). With r c /a < 10~ 2 , equation (4) yields a lower limit for the optical 
depth. The cloud luminosity is then Ld/T* ~ rr 2 /4a 2 . Although r depends on r c , the 
luminosity is independent of cloud geometry. 


Clouds viewed in the orbital plane of the disk can eclipse the central star. The visual 
extinction through the cloud is Ay = 1.086r. The maximum eclipse depth is 

Am ~M/) 2 (£) • (11) 

where i?* is the radius of the central star. At 1 AU, a debris cloud with r c /a ~ 10” 2 in orbit 
around a 1-3 M 0 star produces eclipses with Am ~ 0.1-1.0 A v . 

Figure 1 shows derived luminosities for debris clouds at 1 AU (rocky objects, solid 
lines), 10 AU (rocky objects, dot-dashed lines), and 100 AU (icy objects, dashed lines). 



Larger projectiles produce more luminous debris clouds. In some cases, the largest projec¬ 
tiles completely shatter the target, leading to a saturation in debris production. Complete 
disruption leads to a maximum in the relative luminosity for a given target, illustrated by 
the plateaus in the Figure. This maximum luminosity is roughly 


L'd,max 

u 


3 x 10- 2 

(lirkur ’ 5 


a < 1.6 AU 
a > 1.6 AU 


( 12 ) 


The luminosity of a debris cloud is sensitive to the distance from the parent star. Our 
results suggest a roughly 7 order of magnitude decline in L d /L+ from 1 ALI to 100 AU. This 
large change is a result of the changing angular scale and amount of debris with semimajor 
axis. For a cloud of fixed mass, the angular size of an opaque cloud as viewed from the central 
star varies as a~ 2 . For a fixed eccentricity, the impact kinetic energy scales as a -1 ; thus, 
the debris mass produced from a collision of two bodies varies as a -1 . Because collisions at 
larger stellar distance produce less debris, the mass of the largest remnant increases roughly 
with a. Taken together, these simple scaling laws suggest L d /L± oc a -4 . Complete disruption 
becomes more important close to the central star, reducing the scaling to L d / L* oc a -3 ' 5 

The luminosity of the debris is also sensitive to the bulk properties of the colliding 
bodies (Figure 2). For collisions between the largest objects at 1-10 AU, the Q d relation 
derived from asteroids in our solar system yields a factor of ten smaller luminosity as the 
‘standard’ relation for basalts from Benz & Asphaug (1999). The debris luminosity from 
smaller collisions is unchanged from the standard relation. For the Durda, Greenberg, & 
Jedicke (1998) Q d relation, small changes, ± 0.1, in /3 g yield factor of 3 variations in the 
debris luminosity. At 100 AU, the Q d relation for rubble piles results in much more luminous 
debris clouds for collisions between 100-1000 km objects. However, the luminosity is still 
small, L d / L* < 10 -6 , compared to typical luminosities, ~ 10 -5 to 10~ 3 , of known debris 
disks. 

Close to the central star, debris clouds can produce modest eclipses (Figure 3). At 1 AU, 
collisions between 300 km and larger objects have Ay ~ 0.01-0.1 mag. Although smaller 
collisions produce observable debris clouds in scattered or thermal emission, they produce 
tiny eclipses with Am < 10~ 3 mag. Beyond ~ 3 ALI, eclipses from all debris clouds are 
unobservable with current techniques. 

Our estimates of the collision times indicate that individual collisions between large 
objects are observable (see also Stern 1994; Grogan et al. 2001; Zhang & Sigurdsson 2003). 
At 1 AU, the collision rates range from 1-10 Myr^ 1 for collisions between two 1000 km objects 
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to 10~ 2 yr -1 for 100 km + 100 km collisions. Smaller collisions occur more frequently but 
do not produce significant eclipses or debris luminosity. At 10-100 AU, observable collisions 
between 300 km and larger objects occur every 10 Myr to 1 Gyr. Because the debris cloud 
merges with the disk on 10 4 to 10 5 yr timescales, the chances of observing these debris clouds 
are small. 


2.3. Uncertainties 

The main uncertainties in our estimates are collision physics and the behavior of the 
debris cloud following the collision (e.g., Durda et ah 2004). Because collisions between small 
objects are frequent, the size distribution of small bodies probably is a small uncertainty. 
For optically thin clouds, the cloud geometry sets the optical depth of the cloud but not the 
luminosity. 

For rocky objects, the collision parameters have been well-studied with numerical sim¬ 
ulations (e.g., Benz & Asphaug 1999; Michel et al 2001; Canup 2004a) and observations 
of the debris of large collisions among asteroids (e.g., Durda, Greenberg, & Jedicke 1998). 
The two approaches yield similar optical depths and cloud luminosities for r < 100 km, 
but differ by a factor of ~ 10 for objects with r ~ 1000 km. Because the number of large 
asteroids with r > 500 km is small, the uncertainty in the observational estimates are large 
(Durda, Greenberg, & Jedicke 1998). Observations of sufficient numbers of collisions in 
debris disks may reduce the uncertainty (see below). 

For icy objects, the uncertainties in the bulk properties are much larger. Numerical 
simulations yield Qb ~ 10 6 erg g~ 4 for crystalline ice or snow (Benz & Asphaug 1999); 
simulations of rubble piles suggest Qb < 10 erg g -1 (Leinhardt & Richardson 2002). The 
rubble pile results agree with analyses of simulations of the break up of comets, such as 
Shoemaker-Levy 9 (e.g., Asphaug & Benz 1996) and C/2000 C6 (Uzzo et al. 2001). The 
apparent break in the size distribution of Kuiper belt objects in the outer solar system 
(Bernstein et al. 2004) also favors Qb < 10 3 erg g -1 (Pan & Sari 2005; Kenyon & Bromley 
2004c). 

The size distribution of Kuiper belt objects also favors a small exponent for the gravity 
component of the disruption energy, (3 g < 1.2-1.3 instead of the (3 g ~ 2 expected from simple 
theory (Kenyon & Bromley 2004c). However, current observations cannot distinguish be¬ 
tween the two relations used in Figures 1-3. Improved constraints from the size distribution 
of Kuiper belt objects may yield tests of these relations (Bernstein et al. 2004). 

The evolution of the debris cloud following the collision is an important uncertainty in 
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our estimates. For a cloud with a radial extent r c , the number of secondary collisions with 
material in the background debris disk ranges from less than one per orbit for the largest 
remnant of the collision to a hundred or more collisions per orbit for debris particles with 
radii of 1 km or smaller. Because the debris cloud expands azimuthally at a rate roughly 
comparable to the collision rate, the optical depth through the cloud probably cannot increase 
with time. However, the mass and luminosity may increase. Depending on the surface 
density in the disk, we estimate that secondary collisions can increase the luminosity of the 
debris cloud by a factor of 2-10. This uncertainty has little impact on our main results. If 
the debris cloud luminosity is a factor of 10 larger than our analytic estimates, individual 
collisions between objects with radii of 100 km or smaller are not detectable at 1 AU. At 
10-100 AU, individual collisions between objects with radii of 1000 km and smaller remain 
invisible. 


3. Coagulation Calculations 
3.1. Methods 

To explore the visibility of individual collisions in more detail, we consider numerical 
calculations of the evolution of particle disks around 1-3 M 0 stars. Our calculations for the 
terrestrial zone around solar-type stars show that individual collisions between 100 km and 
larger objects are visible as short-lived increases in the disk luminosity (Kenyon & Bromley 
2004b, see also Grogan et al. 2001). Beyond 30 AU, individual collisions are not visible as 
discrete events in our calculations (Kenyon & Bromley 2004a). Here, we consider a broader 
suite of calculations designed to test the visibility of individual collisions at 1-100 AU. 

To evolve a particle disk in time, we follow the Safronov (1969) statistical approach 
to calculate the collisional evolution of an ensemble of planetesimals in orbit around a star 
of mass M* and luminosity L* (Kenyon & Bromley 2004c, and references therein). The 
model grid contains N concentric annuli with widths 5ai centered at heliocentric distances 
. Calculations begin with a differential mass distribution n(m^) of bodies with orbital 
eccentricity e^t) and inclination i ifc (t). 

To evolve the mass and velocity distributions in time, we solve the coagulation and 
Fokker-Planck equations for bodies undergoing inelastic collisions, drag forces, and long- 
range gravitational interactions (Kenyon & Bromley 2002). We adopt collision rates from 
kinetic theory - the particle-in-a-box method - and use an energy-scaling algorithm to assign 
collision outcomes (Kenyon & Bromley 2004c). We derive changes in orbital parameters 
from gas drag, dynamical friction, and viscous stirring (Adachi et al. 1976; Ohtsuki, Stewart, 
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& Ida 2002). Dynamical friction transfers kinetic energy from large bodies to small bodies 
and drives a system to energy equipartition. Viscous stirring transfers angular momentum 
between bodies and increases the velocities of all planetesimals. 

To compute dust masses and luminosities, we use the dust production rates and particle 
scale heights derived from the coagulation code as input to a dust evolution code (Kenyon 
& Bromley 2004a). This calculation includes collisions between dust grains and Poynting- 
Robertson drag. We use a simple radiative transfer method to derive the optical depth 
through the disk; dust luminosities follow from the scale height derived from the coagulation 
code. 

Our numerical calculations begin with 1-1000 nr planetesimals in 32-64 concentric annuli 
with initial surface density E = E 0 (a/1 AU)” 3 / 2 (see, for example, Kenyon & Luu 1999; 
Kenyon & Bromley 2004a,b). Table 1 lists a in and a out , the inner and outer radii of the 
computational grid, p, the mass density of a planetesimal, E 0 , and t g , the timescale for 
exponential gas removal (Kenyon & Luu 1999). We divide the initial continuous mass 
distribution of planetesimals into a differential mass distribution with discrete mass batches. 
The spacing between successive mass batches is 5 = We add mass batches as 

planetesimals grow in mass. Each planetesimal batch begins in a nearly circular orbit with 
eccentricity eo and inclination i 0 = eo/2. For these calculations, 5 = 1.4-2 and eo = 10” 5 - 
10 ” 4 . 

In the next two sections, we consider the evolution of icy and rocky planetesimals in orbit 
around solar-type and A-type stars. Because our calculations do not include gas accretion, 
we avoid the region of the gas giant planets and restrict our attention to the terrestrial zone 
and the region of the Kuiper belt around solar-type stars. For A-type stars, we consider two 
representative regions at 3-20 AU and at 30-150 ALL Although gas giants may form around 
A-type stars, our calculations provide reasonable estimates for the formation of the rocky 
or icy cores of gas giants. As a planetary core accretes gas from the disk, it will continue 
to stir up the leftover planetesimals and remove the debris more rapidly than estimated in 
our calculations. Thus our estimates for the lifetimes and dust luminosity of debris disks are 
probably reasonable upper limits for these systems. 


3.2. Solar-type stars: results at 0.4—2 AU and 30-150 AU 

At 1 AU, the growth of planetesimals into planets follows a standard pattern (Wetherill 
& Stewart 1993; Weidenschilling et al. 1997). All bodies begin with small orbital eccentric¬ 
ities. Gas drag maintains small e for the smallest bodies; dynamical friction maintains small 
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e for the largest bodies. Thus, collisions produce larger bodies instead of smaller fragments. 
After ~ 10 4 yr, 1000 km and larger bodies form and stir up the orbits of planetesimals 
with radii of 0.1-10 km. Collisions between the leftover planetesimals then produce smaller 
fragments. The largest bodies sweep up some of the fragments, but most of the fragments 
collide with other fragments to produce still smaller fragments. Because the largest bodies 
are nearly immune to fragmentation, this collisional cascade concentrates most of the initial 
mass in the largest and the smallest bodies. 

During the earliest stages of our calculations, the collisional cascade produces an ob¬ 
servable amount of dust (Figure 4). After ~ 10 3 yr, the dust produces an order of magnitude 
more radiation at mid-infrared (rnid-IR) wavelengths than the central solar-type star. Con¬ 
tinued dust production maintains a roughly constant rnid-IR excess for ~ 10 6 yr. Because 
the large bodies are slowly depleted, the dust production rate then slows and the mid-IR 
excess declines. In these calculations, the mid-IR excess becomes comparable to the stellar 
radiation after ~ a few x 10 7 yr. 

In addition to the slow evolution of the ‘background’ mid-IR excess, individual collisions 
produce large spikes in the dust luminosity (Figure 4). During the early stages of the 
evolution, collisions between large objects tend to produce even larger objects and are also 
too frequent for any single, disruptive collision to modify the mid-IR excess significantly. As 
the average dust production rate declines, individual collisions between large objects are less 
frequent but produce more dust. Thus, after ~ 10 6 yr, spikes from large collisions dominate 
the evolution of the mid-IR luminosity. 

To learn more about the nature of these spikes, we consider the dust production rate 
as a function of time (Figure 5). At this point in the evolution (~ 5 x 10 6 yr), numerous 
collisions between small objects set the background dust production rate. Because we use a 
random number generator to assign collision outcomes (Wetherill & Stewart 1993; Kenyon 
& Luu 1998), an occasional lack of collisions leads to a smaller than average dust production 
rate. In addition, occasional collisions between larger objects produce distinct spikes in the 
dust production rate. The frequency distribution of the spikes in the dust production rate 
tracks the size distribution of large bodies; small spikes are much more frequent than large 
spikes. The smaller spikes, where the dust production rate is a few times the background 
rate, produce little change in the mid-IR excess. Larger spikes - produced by collisions 
between 100 km or larger objects - yield fluctuations in the mid-IR excess ranging from a 
few percent to more than 100%. The largest fluctuations in Figure 5 produce the largest 
spikes in the evolution of the mid-IR excess shown in Figure 4. 

At 30 AU, the evolution follows a different path. Because the orbital periods are longer 
and the surface density smaller, planets take much longer to form at 30 ALT than at 1 ALL 



Large bodies begin to form at ~ 10 Myr and soon contain a few percent of the initial 
mass. These objects stir up leftover planetesimals along their orbits and begin the collisional 
cascade. In these models, large bodies grow more slowly, and the leftovers are more easily 
shattered. Collisions then efficiently convert most of the initial mass into smaller objects 
which are removed by radiation pressure and Poynting-Robertson drag. Thus, the final 
ensemble of large bodies contains less than 5%, and sometimes less than 1%, of the initial 
mass in solid objects. 

The evolution of the IR excess is also different at 30 ALI than at 1 ALI (Figure 6). 
Because the dust is farther away from the central star, the dust is cooler, ~ 30-40 K instead 
of ~ 300-400 K. At the start of our calculations, this dust produces a significant IR excess, 
which slowly declines as small bodies merge into larger bodies. Once the collisional cascade 
begins at ~ 10 Myr, the far-IR excess rises significantly. At shorter wavelengths, the dust 
produces a small excess. As in calculations at 1 AU, the cascade maintains a roughly constant 
excess until planet formation reaches the outer edge of the disk at 100-150 AU. Poynting- 
Robertson drag and radiation pressure then remove small grains faster than collisions can 
replenish them. The IR excess then rapidly declines. However, individual collisions never 
contribute a significant amount to the IR excess. In our calculations, the evolution of the 
IR excess is smooth, with l%-2% fluctuations about the general trend. 

In our calculations at 30 AU, there are several reasons for the lack of significant spikes 
in the IR excess. Because the collisional cascade is efficient at removing 1-10 km bodies from 
the evolution, few of these objects grow into 100 km and larger objects. Collisions between 
individual objects also do not produce an observable amount of dust (Figs. 1-2). To test 
the relative importance of these processes, we calculated several models with more robust 
planetesimals, where the collisional cascade leaves more material in larger bodies. These 
calculations suggest that individual collisions never produce observable amounts of dust. 

Our suite of calculations confirms the basic conclusions of our analytic model. At 30 
AU, binary collisions never produce an observable IR excess. Although the dust in our 
calculations produces a large far-IR excess, the evolution of the excess is very smooth. At 
1 AU, occasional collisions between large objects produce observable changes in the mid-IR 
excess. Small, few percent, changes in the mid-IR excess should recur on timescales of tens 
to hundreds of years. Large changes are much less frequent and recur on thousand year or 
longer timescales. 
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3.3. A-type stars: results at 3—20 AU and 30—150 AU 

Calculations for planet formation around A-type stars follow the same pattern as for 
solar-type stars. During the early stages, runaway growth leads to the formation of proto¬ 
planets with radii of 1000 km and larger. At 3 AU, it takes about 10 5 yr for 1000-2000 km 
objects to form. Due to the large gradient in surface density, the growth of large bodies 
slowly propagates to larger distances (Lissauer 1987). The timescale for the formation of 
1000 km objects is ~ 10 6 yr at 9 AU and ~ 10' yr at 18-20 AU. 

Once large planets form, their gravity stirs up the leftover planetesimals along their 
orbits. Collisions then produce fragments instead of mergers, and the collisional cascade 
begins. In simulations at 3-20 AU, the dust luminosity starts to rise at ~ 10 4 yr and 
peaks at ~ 10 6 yr (Figure 7). As planet formation moves from ~ 10 AU to ~ 20 AU, the 
dust produced in the collisional cascade has smaller optical depth and subtends a smaller 
solid angle as seen from the central star. Thus, the dust luminosity declines. Once planet 
formation reaches the outer edge of the 20 AU disk, Poynting-Robertson drag removes dust 
more rapidly than collisions can replenish it. The dust luminosity then declines more rapidly 
and falls to very small levels in 10 9 yr. 

In calculations at 30-150 AU (Kenyon & Bromley 2004a), planet formation leads to 
dust production on much longer timescales. In these models, large planets begin to form at 
the inner edge of the disk at 10-20 Myr. After 1 Gyr, planet formation reaches the outer 
edge of the disk. Throughout this period, the collisional cascade maintains a large reservoir 
of small dust grains, which produce a measurable dust luminosity (Figure 7). From 10 Myr 
to 1 Gyr, the dust luminosity smoothly declines from L d /L 0 ~ 10 3 to L d /L 0 ~ 10 4 . After 
1 Gyr, Poynting-Robertson drag rapidly removes dust and the dust luminosity declines more 
rapidly. 

Figure 8 shows the evolution of the 25 //m excess for a planetesimal disk around an 
A-type star with mass M* = 3 M 0 , luminosity L* = 50 L 0 , and temperature T* = 9500 K. 
In both cases, the excess first declines with time as small bodies merge into larger and larger 
bodies. Once the collisional cascade begins, the 24 /mi excess rises. Small fluctuations from 
individual collisions mark the rise of the excess at 3 AU. At 30 AU, the rise is smooth. At the 
peak of both calculations, small fluctuations in the excess are apparent. At 3 AU, individual 
collisions produce the small fluctuations. At 30 AU, the occasional merger of large objects 
produces more stirring and a temporary increase in the dust production rate. 

These results also confirm the basic predictions of the analytic collision model. In 
planetesimal disks surrounding A-type stars, binary collisions at 3-5 AU produce modest, 
~ 10%, fluctuations in the mid-IR excess. Farther out in the disk, binary collisions do not 



produce enough dust to change the amount of mid-IR excess. However, mergers of larger 
objects can produce modest fluctuations in the mid-IR excess by increasing the local stirring 
rate, which leads to more dust and a larger mid-IR excess. 


4. N-BODY CALCULATIONS 

When binary collisions occur, our model predicts two observable changes in the lumi¬ 
nosity. After a collision produces a debris cloud, we expect a rapid rise in infrared brightness. 
As orbital shear spreads the debris cloud into a ring, secondary collisions with material in 
the rest of the disk remove smaller particles. The mid-infrared flux then declines. In our 
picture, this decline occurs on the timescale for the cloud to spread into a ring, ~ 10 2 orbital 
periods. 

For edge-on systems, we expect the debris cloud to eclipse the central star. The predicted 
eclipse depth ranges from ~ 0.5 mag for large collisions to less than 0.01 mag for smaller 
collisions. As the debris cloud spreads into a ring, the eclipses should last longer and become 
shallower. This behavior distinguishes a debris cloud from the transit of an exosolar planet, 
where the eclipse depth and duration are constant (see, for example, Alonso et al. 2004, 
and references therein). 

To illustrate the appearance of eclipses from a debris cloud, we use IV-body simulations 
to follow the evolution of debris following a major two body collision. For these examples, 
we consider collisions between equal mass bodies with radius 300 km and 1000 km. 

To begin the evolution, we place 10,000 “dust” particles in a small volume of radius 
r c /a = 10 -3 around a guiding center on a circular orbit around a star of mass M* = 1 
Mq. We give each dust particle a relative velocity V/Vk ~ 10~ 3 . In one configuration, the 
positions and velocities have Gaussian random distributions (generating a “random blob”). 
In another configuration, the particle density declines as r~ 2 out to radius r c , with velocities 
proportional to the distance from the guiding center (yielding a homologous “explosion”). 
These two configurations represent two reasonable choices for the initial velocity distribution 
(e.g., Michel et al 2001). 

Each IV-body particle represents a population of dust with mass distributed according 
to equation (10), with a = 3.5. The minimum size of the dust particles is 1 /jrn. We set 
the maximum size in one of two ways. In the first case, we assume that all of the mass in 
the collision goes into “dust particles” with a single-index power law size distribution and 
a single large remnant containing 10%-15% of the mass. In the second model we suppose 
that 20% of the mass ends up in a single large object, with the remaining 80% in dust 
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with a maximum radius of 10 km. Although these choices do not exhaust the possible mass 
configurations following a completely disruptive collision, they provide a good illustration of 
the range of eclipse depths and evolution with time. 

To derive the properties of the eclipses, we calculate the position of the particles on 
Keplerian orbits evaluated with an 8 th -order ODE solver. We divide the surface of a star 
with radius 1 R e , into 0(100) elements and project the particle positions onto the stellar 
disk. For each path from the observer through the particles to the stellar disk, we compute 
the optical depth r(t) and the emergent intensity I'(t) = e~ T ^I^ where the unobscured 
intensity from the star is / 0 = I t . This procedure assumes that the mass distribution for 
all particles is identical. We then derive a light curve from -2.5 log QT) I[{t))/I q. 

Figure 9 illustrates light curves for the collision of two 1000 km objects. The lower 
(upper) panel shows results for a random blob (explosion) when a large remnant contains 
20% of the initial mass. In both simulations, the eclipse depth is ~ 0.5-0.7 mag ~ 6 months 
after the collision. Each eclipse lasts ~ 5-10 days, ~ 2% to 4% of the orbital period. After 
10-20 orbits, the eclipses are ~ 10% of their original depth and last ~ 10 times longer than 
eclipses immediately after the collision. After ~ 100 orbits, dust from both the blob and 
the explosion lies in a roughly uniform ring around the star. Although tiny eclipses remain 
visible, most of the dust produces a uniform extinction of ~ 2% to 3% of the stellar flux. 

In these simulations, dust from the random blob expands more rapidly and merges into 
the background debris disk faster than dust from the explosion. This behavior is a feature of 
the initial conditions. In the homologous explosion, the dust particles are less concentrated to 
the center, which produces shallower eclipses relative to the random blob. However, particles 
close to the center of the explosion expand more slowly than those in the blob. Thus, these 
particles remain closer to the point of the collision and continue to produce narrow, shallow 
eclipses when most of the rest of the blob has expanded into a ring. 

Figure 10 repeats Figure 9 for the collision of two 300 km bodies. Because these collisions 
produce less dust, the eclipses have roughly the same duration but are shallower. As for the 
collision of 1000 km bodies, the eclipses slowly fade away and last a longer and longer fraction 
of the orbital period. After 30-100 orbits, orbital shear transforms the debris cloud into a 
ring, which obscures ~ 1% to 2% of the stellar surface. 

Figure 11 compares results for configurations where we place all of the debris in “dust” 
particles with radii of 10 km or smaller. The lower (upper) panel shows light curves for the 
collision of two 1000 km (300 km) bodies. Compared to configurations where 20% of the 
initial mass in contained in a single large object (Figures 9-10), the optical depth in the 
debris cloud is larger and eclipses are much deeper. The timescale for the cloud to expand 
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into a ring is identical to the timescale derived for the other cases. Thus, the magnitude of 
the eclipses declines over 30-100 orbital periods. 

These calculations illustrate the simplest evolution of a debris cloud formed in a large 
binary collision. Calculations with more complex initial configurations produce a remarkable 
variety of eclipse light curves. Shepherd moons or planets can confine the debris cloud and 
thus can slow the decline of the eclipse depth with time. Planets in orbits inclined to the 
orbital plane of the debris cloud lead to a precessing cloud or ring with very different eclipse 
morphology. Winn et al. (2004) and Chiang & Murray-Clay (2004) show how the eclipse 
morphology changes when an opaque cloud or ring orbits an eccentric binary star and apply 
their results to the young variable star KH 15D (Hamilton et al. 2001). We derive similar 
results when a debris cloud spreads into a ring. Because the range of possible geometries 
is large, the discovery of additional eclipsing systems with rings or blobs of dust will guide 
future numerical calculations. 


5. DISCUSSION AND SUMMARY 

During the late stages of planet formation, binary collisions between leftover, large 
planetesimals in orbit around solar-type stars are potentially observable. At 1 AU, these 
collisions produce large spikes, ~ 1 mag or less, in the mid-infrared luminosity (Figure 4). In 
systems viewed edge-on, the debris cloud eclipses the parent star (Figures 9-11). Maximum 
eclipse depths range from ~ 0.5 mag for collisions between 1000 km objects to ~ 0.01 mag 
or less for collisions between 100-300 km objects. Typical eclipses last 5-10 days. As the 
debris cloud expands into a ring, the eclipses fade and grow longer. The lifetime of the debris 
cloud as a distinct entity is ~ 10 2 orbital periods. 

At large distances from the parent star, > 3 AU, binary collisions are more difficult 
to detect. Because collisions are less frequent and produce less dust, they are harder to 
detect against the emission from the background debris disk. Our results indicate that the 
maximum luminosity from individual collisions declines as a ~ 3 5 (equation 12). Wyatt & 
Dent (2002) derive similar conclusions from different approaches. 

To estimate the likelihood of detecting binary collisions, we consider Keplei 3 , a satel¬ 
lite designed to detect the transits of exosolar planets across their parent star. During a 4 
yr mission, the Kepler team plans to acquire very high precision photometry for 10 5 stars. 
Roughly 0.5% of these stars have favorable orbital inclinations to detect transits from terres- 


3 http: //www.kepler.arc.nasa.gov 
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trial planets. For the planned observing frequency of once every few hours and an expected 
precision of 0.01% or better, the expected detection rate is ~ 50-100 terrestrial planets per 
10 5 stars. 

To modify this estimate for the detection of debris clouds, we use the results from our 
coagulation and n-body calculations. Because the scale height of the debris disk is larger than 
the radius of a terrestrial planet, ~ 2% to 3% of potential solar systems have orientations 
that produce eclipses. If the lifetime of the collision phase is a few hundred Myr, ~ 10% of 
Kepler targets are young enough to observe collisions. Thus, eclipses from debris clouds are 
possible in ~ a few hundred stars of the Kepler sample. 

To estimate the probability of detecting binary collisions in these stars, we adopt the 
collision rates from the coagulation calculations and the lifetimes derived from the n-body 
results. For the anticipated precision of Kepler, binary collisions between 100 km objects 
yield a detectable eclipse. With a collision rate of 10~ 2 yr' 1 for these objects, the probability 
of a collision during the Kepler lifetime is roughly 100%. In a 4 yr mission, the eclipse depth 
can change by a factor of two or more. The duration of each eclipse lengthens to maintain 
a roughly constant total optical depth, integrated over the complete eclipse. 

The prospect for observing collisions between the largest objects is poorer but still 
significant. The rate for collisions between two 1000 km objects is < 10~ 5 yr' 1 . This 
rate yields a ~ 1% probability for detecting one large collision during the Kepler mission. 
However, with a typical lifetime of 30-100 yr, the probability of observing eclipses after a 
large collision is ~ 10 times larger than observing a collision during the mission. Thus, it 
is possible that Kepler will detect eclipses from the aftermath of one large collision between 
two 1000 km objects. These eclipses should lengthen and decrease in depth throughout the 
4 yr Kepler mission. 

Gronnd-based surveys could also detect collisions in the terrestrial zone. The OGLE 
and TRES projects have identified several transiting exosolar planets (Alonso et al. 2004; 
Konacki et al. 2003; Pont et al. 2004). The PASS project envisions a complete survey 
of planetary transits for all stars with V ~ 5.5-10.5 (Deeg et al. 2004). As part of a 
much deeper all-sky survey, Pan-STARRS 4 might also detect eclipses from debris clouds. To 
detect the unique signal from a large two body collision, these projects need to distinguish 
real events with diminishing depth and lengthening duration from spurious signals and from 
repetitive eclipses with constant depth and duration. 

Successful detections of the dust produced in large two body collisions can teach us about 


4 http: //pan-starrs.ifa.hawaii.edu/public/index.html 
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terrestrial planet formation. Characterizing the frequency of events associated with stars of 
different ages and spectral types yields estimates for the space density of large bodies as a 
function of stellar mass and age. Simultaneous optical, infrared, and subnnn observations 
provide measures of the dust albedo and size distribution. Coupled with direct detection of 
exosolar terrestrial planets, measurements of the properties of the terrestrial environment 
allow a better understanding of the formation of the Earth-Moon system and, ultimately, 
the formation of the current environment on the Earth and other terrestrial planets. 
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Fig. 1.— Luminosity of debris clouds as a function of projectile radius at 1 AU (solid curves), 
10 AU (dot-dashed curves), and 100 AU (dashed curves). The log of the target radius (in 
cm) is listed to the right of each solid curve. The target radii for the other curves follow 
the same pattern. The curves assume bulk properties of (a) basalts from Benz & Asphaug 
(1999) at 1-10 AU, and (b) modified ices from Benz & Asphaug (1999) as listed in the 
main text. At 1 AU, collisions between 100 km and larger objects produce observable debris 
clouds. At 3 AU and beyond, individual debris clouds are not visible. 
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Fig. 2.— As in Figure 1 for objects with (a) bulk properties for asteroids from Durda, Green¬ 
berg, & Jedicke (1998, solid curves), and (b) bulk properties of rubble piles from Leinhardt 
& Richardson (2002, dashed curves). 
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Fig. 3.— Optical depths for the models in Figure 1 
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Fig. 4.— Evolution of the 24 /jrn excess as a function of time for planetesimal disks at 1.0 
AU. The radial extent of the disk is 0.68-1.32 AU (0.4-2.0 AU) for the dashed (solid) curve. 
In both cases, formation of lunar mass objects leads to the production of copious amounts 
of dust and a significant excess of radiation at 24 //m >1)- As objects grow 

into Mars-sized or larger planets, less frequent collisions produce smaller amounts of dust. 
Individual large collisions are then responsible for large increases in the dust excess. 
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Fig. 5.— Variation of the dust production rate with time for the 0.4-2 AU calculation in 
Figure 4 at 5 x 10 6 yr. Numerous collisions between 1 m to 1 km objects produce the roughly 
constant dust production rate of ~ 5 x 10 18 g yr^ 1 . Binary collisions produce occasional 
spikes in the dust production rate. 
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- Evolution of the IR excess at 60 /mi and at 450 /mi for a planetes: 
VU around a solar-type star. 
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Fig. 7.— Evolution of the dust luminosity for models at 3-20 AU and at 30-150 AU. The 
central A-type star has M* = 3M 0 , L+ = 5OL 0 , and T* = 9500 K. Binary collisions produce 
the small variations in dust luminosity at the peak luminosity of each model. During the 
decline, individual collisions are not visible as distinct luminosity spikes. 
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Fig. 9.— Variation in brightness of a solar-type star eclipsed by a debris cloud formed from 
the collision of two 1000 km objects at 1 AU. In these calculations, a large remnant contains 
~ 15% of the initial mass; the remaining mass is in smaller objects with a power-law size 
distribution (q = 3.5). The legend in each panel indicates the time after the collision for 
each curve. Orbital shear and internal expansion of the cloud produce the gradual evolution 
in eclipse depth and duration. 
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Fig. 11.— As in Figures 9-10 for collisions where 20% of the mass is in a single large remnant 
and 80% of the mass is in objects with radii of 10 km or smaller. 






Table 1. Model Parameters 


Parameter 

Symbol 

G Star Models 

A Star Models 



A 

B 

A 

B 

Stellar mass (M 0 ) 

M* 

1 

1 

3 

3 

Stellar luminosity (L 0 ) 

u 

1 

1 

75 

75 

Inner radius (AU) 


0.4 

30 

3 

30 

Outer radius (AU) 

&out 

2.0 

150 

20 

150 

Particle mass density (g cm -3 ) 

p 

3.0 

1.5 

3.0 

1.5 

Surface density at 1 AU (g cm -2 ) 

So 

7.5-30 

15-60 

5-50 

5-50 

Gas removal timescale (yr) 

tg 

10 6 -10 7 

10 7 

10 6 -10 7 

10 7 






